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Abstract. In this paper, we demonstrate that the modified Emden type equation 
(MEE), x + axx + (3x 3 — 0, is integrable either explicitly or by quadrature for any 
value of a and (3. We also prove that the MEE possesses appropriate time-independent 
Hamiltonian function for the full range of parameters a and f3. In addition, we show 
that the MEE is intimately connected with two well known nonlinear models, namely 
the force-free Duffing type oscillator equation and the two dimensional Lotka-Volterra 
(LV) equation and thus the complete integrability of the latter two models can also be 
understood in terms of the MEE. 



PACS numbers: 02.30.Hq, 02.30.1k, 05.45.-a 
1. Introduction 

One of the well discussed models in nonlinear dynamics is the modified Emden type 
equation (MEE), also called the modified Painleve-Ince equation, 

x + axx + (3x 3 = 0, (1) 

where over dot denotes differentiation with respect to time and a and (3 are arbitrary 
parameters. This equation has received attention from both mathematicians and 
physicists for more than a century [1-5]. For example, in the nineteenth century, 
Painleve had studied this equation and identified general solution for two parametric 
choices, namely, (i) (3 — ^ and (ii) /3 = —a 2 [1,2]. The above differential equation 
(CQ) arises in a variety of mathematical problems such as univalued functions defined by 
second order differential equations [6] and the Riccati equation [7]. On the other hand 
physicists have shown that this equation arises in different contexts: for example, it 
occurs in the study of equilibrium configurations of a spherical gas cloud acting under 
the mutual attraction of its molecules and subject to the laws of thermodynamics [8, 9] 
and in the modelling of the fusion of pellets [10]. It also governs spherically symmetric 
expansion or collapse of a relativistically gravitating mass [11]. This equation can also 
be thought of as a one-dimensional analogue [7,12] of the boson 'gauge-theory' equations 
introduced by Yang and Mills. Apart from the above, for the past two decades or so the 
invariance and integrability properties of this equation alone have been studied in detail 
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by a number of authors, see for example Ref. [13-23]. In a nutshell, the MEE ([T]) has 
been found to possess explicit general solution only for the following parametric choices, 
that is, (i) a = 0, (ii) (3 = 0, (Hi) (3 = ^ and (iv) (3 = —a 2 . While the cases (i) and 
(ii) can be integrated trivially, in the third case the equation is linearizable to a free 
particle equation and in the fourth case the general solution can be expressed in terms 
of Weierstrass elliptic function [1-5,13-24]. Equation ([I]) has also been noted to possess 
Painleve property only for certain values of r = fs(a. ± \J a 2 — 8(3) [17, 19]. Finally, 
we mention that equation (T1J admits a two parameter Lie point symmetry group for 
arbitrary values of a and (3, while for the chioce /3 = %- the equation JT]) possesses eight- 
parameter Lie point symmetries [13]. However, the general solution for the equation (TjQ) 
with a and (3 are arbitrary is yet to be explored. 

Very recently [23], the present authors have studied a generalized version of this 
equation from a different perspective and shown that the equation (jTJ, for a 2 > 8/3, 
possesses time independent integrals and admits a Hamiltonian formalism which in 
turn ensures its complete integrability [23]. However, due to the complicated form of 
the first integral the general solution was not obtained in the previous work. Keeping 
in mind the historical importance and popularity of this model we kept exploring the 
general solution of equation ([!]) in phase-space. Based on our investigation, in this 
paper, we construct the time independent integrals for the equation ([I]) for arbitrary 
values of a and (3 (including the case a 2 < 8(3). Since the first integrals are not in 
simple polynomial forms it is difficult to obtain the general solution by just directly 
integrating them. In order to overcome this difficulty, first we identify time independent 
Hamiltonians from these time independent integrals and making use of suitable canonical 
transformations we convert the Hamiltonians into standard forms. We then integrate 
the new Hamiltonians and obtain the general solutions or reduce to quadratures. In this 
way we report the general solution for the equation ([1]) for arbitrary values of a and (3 for 
the first time. Our motivation to explore this solution is, as we see below, also due to the 
fact that the MEE ([I]) is not a stand-alone model. It is intimately connected to force-free 
Duffing oscillator type equation (vide equation (13T1) below) and two dimensional Lotka- 
Volterra (LV) equation (vide equation fl33|) below). The popularity and importance of 
these two models need no emphasis [23,25-29]. Thus exploring the general solution for 
the equation ([[]) also serves to establish the complete integrability of these two models 
for appropriate parameters besides understanding the dynamics of the other models 
mentioned in the introduction. 

The plan of the paper is as follows. In Section [2], we give the time independent 
integrals and corresponting Hamiltonians for the MEE (OQ). Using suitable canonical 
transformations, we obtain the general solutions for the parameter ranges a 2 = 8(3, a 2 > 
8/3 and a 2 < 8(3 separately in Section [3j In Section HJ we show that the MEE ([1]) is 
intimately connected with two other well known nonlinear models, namely, the force-free 
Duffing oscillator type equation and two dimensional Lotka-Volterra equation. Finally, 
in Section [5] we summarize our results. 
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2. Time independent integrals and Hamiltonian description 

As pointed out in the introduction, the MEE (OQ) cannot be straightforwardly integrated. 
Making use of the modified Prelle-Singer method developed by us recently [22, 23] 
following the earlier work of Duarte et al [30], in this section we first identify the first 
integrals separately for each of the three ranges (i) a 2 = 8(3, (ii) a 2 > 8(3 and (iii) 
a 2 < 8(3. Then we identify suitable canonical Hamiltonian description for each of these 
cases. 



2.1. Time independent integrals 

In Ref. [23], the time independent integrals for the equation (QQ) with a 2 > 8(3 has been 
reported using the modified Prelle-Singer procedure. However, improving the ansatz 
given in Ref. [23] one can obtain the time independent integrals for all values of a and 
(3 and the method of deriving the integrals for the equation (CQ) with a and (3 arbitrary 
is given in Appendix A. Using this method we identify the following time independent 
first integrals for the equation (CQ) with a and (3 arbitrary, that is, 

Case 1: a 2 = 8(3 (r = 2) 

/ = log(aV + Aax) - ^ \ , (2) 

a 2 x 2 + Aax 

Case 2: a 2 > 8(3 (r ^ 0,1,2) 

r (r-1)/. (r-1) 2 Y- r {. a 2 \ 

1=) J - [x + - V [x + — x 2 ) 3 

r-2 V 2r \ r V ; 



Case 3: a 2 < 8(3 



- log(2x 2 + ax 2 x + (3x A ) H — tan 1 

2 LU 



ax + 2(3x 



2i 



LUX 



(4) 



where to = a/8/3 — a 2 and r = fz(ot ± \J ot 2 — 8fi). Note that r = and 1 correspond 
to the trivial cases a = and (3 = 0, respectively, and so they are not considered here 
separately. 

As it is very difficult to integrate equations (G3)-(J4]) and obtain the general solutions 
by direct integration, we correlate these integrals with appropriate Hamiltonians. In the 
following we briefly give the method of obtain the Hamiltonian from the known time 
independent integral. 

2.2. Hamiltonian description 

To explore the Hamiltonian description of (j2J), let us assume the existence of a 
Hamiltonian 

I(x, x) = H(x, p) = px — L(x, x), (5) 

where L is the Lagrangian and p is the canonically conjugate momentum. Then 
dl OH dp . dL dp . 
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From 



we identify 

P = 



—dx 

x 



/(*), (7) 

where f(x) is an arbitrary function of x and the Lagrangian L = px — I(x,x). Here, 
without loss of generality, we take f(x) = 0. Substituting the integrals ©-(jl]) into (|7j) 
and integrating the resultant integrals we can obtain the expression for the canonical 
momentum p. Substituting back the latter into the equation (J5D and simplifying the 
resultant equation we arrive at the following Lagrangian, 



L 



and Hamiltonian 



H 



-log(x + fx 2 ) 

1 (± + tlll ax 2)(2-r) 
(2-r) \ ' 2r / ' 



a" = 
a 2 > 



IP 



tan 



Ax+ax 2 



'—) — atan 1 



ax+2l3x 2 



— \ log(2x 2 + ax 2 x + (3x 4 



(8) 



a 2 < 8p 



(r 



-(p) r " 



fpx 2 , 

r-2 (? 
1 — _ 



2r 



1) 2 
— ax p, 



a 



a 2 > 



(9) 



-\og[x 4 sec 2 (^x 2 p)} - % 2 p, 
where the canonically conjugate momentum 



a 2 < A(3 



P 



2r 



1) 2 

-ax 



l-r 



UJX^ 



-tan 



Ax + ax 
2lox 2 



or 



o 2 > 



o 2 < 



(10) 



respectively for the MEE ([T]). 

One can easily check that the canonical equations of motion for the above 
Hamiltonians are nothing but the equation of motion ([1]) of the MEE in the appropriate 
parametric regimes. 

3. General solutions 

In this section, we consider each of the above three cases separately, namely, (i) o 2 = 8/3, 
(ii) a 2 > 8f3 and (iii) o 2 < 8/3, and obtain their respective general solutions using suitable 
canonical transformations. 



3.1. Case 1: a 2 = 8/3 

To derive the general solution for this parametric choice first we consider the 
Hamiltonian given in (jHJ) for o 2 = 8/3 as 



H = log(- 



4a, 
~P' 



a 



-px 



(11) 
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Introducing the canonical transformation 
4P all 2 



x 



aU'' 



P 



(12) 



(13) 



the Hamiltonian (fTTl) can be recast into the standard form 
H=\P 2 + log(§) = E, 

where E is a constant. From the corresponding canonical equations U = P, P = jj, 
equation ( flBl can be rewritten as 



E=±U 2 -2log(U), £ = £-log(32). 



Rewriting the above equation we get 
dU 

-. = dt. 



'2E + 4\og(U) 
Integrating the equation ([TBI) we obtain 

-\[E+2eir 1 (z) 2 
U = e V 



(14) 



(15) 



(16) 



whore : = ^.J" ' - and to is the second arbitrary integration constant and erf is the 



7T 



error function [31]. 

Substituting the equation (TlBT) into equation (Tl2l) we get the general solution for the 
equation ([1]), with the parametric choice (3 = in the form (after some modifications) 



8 _ if^err 1 ^) _ 2ef(t -t) 
— ert ( z)p. \ / . z = = 



x(t) = —erf 1 (z)e 
a 



(17) 



In Fig. 1, we have plotted the solution of MEE (pj) given by the expression (fTTj) for the 
parametric choice a 2 = 8/? for different initial conditions. 
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Figure 1. Solution plot for case 



a 2 = 8/3 for different initial conditions. 
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3.2. Case 2: a 2 > 8(3 

Now let us consider the Hamiltonian for the case a 2 > 8(3 from ([9]), that is, 

H = " (18) 

(r — 2) 2r 

Interestingly, here we identify a canonical transformation for the Hamiltonian (JT8l) 

in the form 

P U 2 

x = -a-, P = -, (19) 
where a = , 2r . It is straightforward to check that when U and P are canonical so do 

a(l— r) ° 

x and p (and vice versa) so that the Hamiltonian H in Eq. (j!8p can be rewritten as 

# = ip 2 - A;[/ m = P, (20) 
2(r-2) 

where m = — — and k = - — -yn-, r- From the canonical equations then we have 

(r - 1) (2a) 2 (r - 2) 

U = P and P = mkU 171 - 1 . 

To obtain the general solution, we introduce a transformation X = U m in equation 

f l20|) so that the latter can be brought to a quadrature of the form 

r x'-^rdX , , . 

'-W wSCTP m ^ ' 2 ' 4 - (21) 

Now fixing = n, the above integral leads us to the following expression [32], that 
is, 



r X n dX _ 1 (X n V2E + 2kX 2En r X n ~ L dX \ 
~ °~ J my/2E + 2kX "ml 2A;(n + |) 2A;(n + |) 7 ^2^ + 2^X7' ( ""' 
On the other hand fixing = — n in (T2T1) . we end up with the following expression [32] 

dX 



t-t, 



o 



mX n ^2E + 2kX 



1 / V2E + 2kX k(n - f ) r dX 



f23) 

m\ 2EX n - 1 (n-l) E(n-l)J X n ~ 1 y / 2E + 2kX!' 

For integer values of n the integrals ( |22i) and (1231) can be integrated explicitly 

by repeated use of these formulas. When n is a noninteger value the final term in 

the integrals ( 1221) and ( 1231) can be integrated in terms of logarithmic function or beta 

function [32] (since the value of n — 1 in the final integral lies between to 1). We also 

note here that for the special choices n = or — | (which corresponds to the case (3 = ^) 

and n = — | or — | (which corresponds to the case (3 = —a 2 ) respectively in equations 

( |22~1) and ( |23l) . one can get the respective solutions reported in the literature [1-5,13-24]. 

3.3. Case 3: a 2 < 8(3 

Finaly, let us focus our attention in the regime a 2 < 8(3. In this parametric regime the 
Hamiltonian takes the form (vide equation (Q) 

H = - log[x sec (—x p)\ — —x p. (24) 
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For the present case, let us choose the canonical transformation in the form 
UP 2 

and transform the Hamiltonian to 

H = \\o g [^sec 2 ^U 2 )]-^U 2 ^E. (26) 

Now making use of the canonical equations, U = — -p and P = ^j((ct — 
uj tan(^[/ 2 ))[/ 2 — 8), equation ([2"B1) can be rewritten as 

E=~ \og[—sec\-U 2 )] - -U 2 . (27) 
Introducing now the transformation V — ^ in fl27|) we arrive at 

^=^og[^ S ec 2 (^)]-^y. (28) 
Integrating the equation (T2"8l we get 

e 2 /• / cj " . — — v 



t-t = —J ^sec(-V)e 8 dV, (29) 
where to is the second integration constant. By substituting W = e v in (1291) we obtain 



where gi = ~"+^ ; q 2 = ^ and F(a, /3; 7; is the hypergeometric function [31,32]. 
4. Connection to two other nonlinear models 

In this section, we show that the MEE ([T]) is intimately connected with two other well 
known nonlinear models, namely, the force-free Duffing oscillator type equation and two 
dimensional Lotka-Volterra (LV) equation. 

4-1. Force-free Duffing oscillator type equation 

The MEE equation (TjQ) can be transformed to the following oscillator equation 

w" + (aw + 7 W + (3w 3 + ^w 2 + = 0, ('= 4~) (31) 

3 9 dr 

through the invertible point transformation x = wei T ,t = — -e~? T , where 7 is an 

arbitrary parameter. Equation fl3T|) includes force-free Duffing oscillator (in the specific 

case a = 0) [23,25,26] and quadratic oscillator with (3 — [23]. 

One can deduce the general solution of the equation f[3"Tj) from the general solution 

2 

of MEE. For example, in the parametric choice j3 = ^- the general solution of the 
equation (13TT) can be derived from (iT7j) in the form 

w(t ) = -erT 1 (2)e V 3 /, (32) 
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where now z = 1= and E and tq are arbitrary integration constants. 

7V 71 " 



Similarly, one can fix the general solution for the equation ( fSTl) in the other parametric 
regimes also, that is, a 2 > 8/3 and a 2 < 8(3. 

4-2. Two dimensional Lotka-Volterra equation 

Let us consider the two dimensional Lotka-Volterra equation of the form 

x = x(ai + a 2 x + a 3 y), y = y(b x + b 2 x + b 3 y), (33) 

where a^s and frj's are six real parameters. Equation (!33l models two species in 
competition in ecology and is being analyzed for the past three decades or so in 
mathematical biology [27-29]. Interestingly, equation (133]) can also be transformed to 
( I3TT) as follows. Rewriting the first equation in ( 1331) for the variable y and substituting 
it into the second equation in equation (!33l we get the following second order ODE for 
x, namely, 

'x — (H ) h ((2a 2 a 2 - 6 2 )x + (2a x h))x 

a 3 x a 3 a 3 

+ {b 2 a 2 — —al)x 3 + (a 2 bi + b 2 a\ — 2a\a 2 —)x 2 + (a±bi — —a\)x = 0. (34) 
a 3 a 3 a 3 

Let us choose the parameters in (1311) in the form b 3 = —a 3 and b\ = ai so that equation 

( |34|) can be brought to the form 

x — ((3a 2 + b 2 )x + 3ai) i; + a 2 (a 2 + b 2 )x 3 + ai(3a 2 + b 2 )x 2 + 2a\x = 0. (35) 

The associated LV equation takes the form 

x = x(ai + a 2 x + a 3 y), y = y(a\ + b 2 x — a 3 y). (36) 

Now comparing the equations (|3~TT) and fl35l) we obtain a = — (3a 2 + b 2 ), (3 = a 2 (a 2 + b 2 ) 
and 7 = —3a!. Choosing a 2 = b 2 {(3 = ^-) the general solution for the LV equation 
( 135]) can be obtained from (!32|) and using this in the first equation in (!36|) we arrive the 
general solution for the LV equation (1361) in the form 

| ( E+2ait+2erf~ 1 (z) 2 



2 

e 



2 _-, | ( £+2ait+2erf 1 (z) 2 

—erf (2)e V . „,. , 

a 2 a 3 ert (zj 



x(t) = —err 1 (z)e 2 \ ' " V, y(t) = . ( 37 ) 



2e^(t(j — e ait ) 

where z = ■= and E and to are arbitrary integration constants. Finally, we 

aiV 71 " 

mention that the general solution for the equation (!35|) for other parametric choices 
can be obtained from equation (I3T1) in a similar manner. Again (to our knowledge) the 
complete integrability of the LV system (1361) in these regimes is new to the literature. 



5. Conclusions 



In this paper, we have shown that the MEE ([T]) is integrable either explicitly or by 
quadratures for any value of a and (3. We have also obtained the time independent 



General solution for the modified Emden type equation 



9 



Hamiltonians for the equation ([T]). We have transformed the Hamiltonians into simpler 
forms, with appropriate canonical transformations, and deduced the general solution by 
direct integration so that our approach helps to understand the dynamics of equation (TjQ) 
in phase-space clearly. Further, we have demonstrated that the complete integrability of 
equation (CQ) also helps one to understand the dynamics of two other nonlinear models, 
namely, the generalized oscillator equation and the two dimensional LV equation. As 
a consequence the solutions which we have explored for the equation ([1]) also provide 
solutions for these models as well. 
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Appendix A. Method of deriving integrals of motions 

In the following we briefly explain the generalized extended or modified Prelle-Singer 
(PS) procedure for second order ODEs [22,23,30] which is used to identify the integrals 
of motions ©-flU). 

To begin with, let us rewrite equation ([I]) in the form 

x = — (axx + /3x 3 ) = (p(x, x). (A.l) 

Further, we assume that the ODE (lA.lj) admits a first integral I(t, x, x) = C, with C 
constant on the solutions, so that the total differential becomes 

dl = I t dt + I x dx + I x dx = 0, (A. 2) 

where each subscript denotes partial differentiation with respect to that variable. 
Rewriting equation (lA.lj) in the form (fidt — dx = and adding a null term S(t,x,x)x 



dt — S(t, x, x)dx to the latter, we obtain that on the solutions the 1-form 

(cf) + Sx)dt — Sdx — dx = 0, cj) = — (axx + (3x 3 ). (A. 3) 

Hence, on the solutions, the 1-forms (1A.2I) and (1A.3|) must be proportional. Multiplying 
f TA~3l by the factor R(t, x, x) which acts as the integrating factors for (1A.3I) . we have on 
the solutions that 

dl = R((f) + Sx)dt - RSdx - Rdx = 0. (A.4) 

Comparing Eq. (1A.2I) with (1A.4I) we have, on the solutions, the relations 

I t = R((j) + ±S), I x = -RS, I x = -R. (A.5) 

Then the compatibility conditions, I tx = I xt , It± = I±t, Ixx = I±x, between the Eqs. ( 1A.5j) . 
provide us 

S t + xS x + 4>S x = - 4> x + 4> X S + S 2 , (A. 6) 
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Rt + xR x + <f>R x = - {fe + S)R, (A.7) 
R x — SR± — RS± = 0. (A. 8) 

Solving equations flA.6j) - f)A.8l) one can obtain expressions for S and R. It may be 
noted that any set of special solutions (S, R) is sufficient for our purpose. Once these 
forms are determined the integral of motion I(t, X j X^j Cclll be deduced from the relation 



n - r 2 



R+ Tx {n - r2) 



where 



Tl 



dx, (A. 9) 



J R((j) + xS)dt, r 2 = J '{RS + -^n)dx 



Equation ( 1A.9I) can be derived straightforwardly by integrating the equation ( 1A.5I) . 

As our motivation is to explore time independent integral of motion for the 
equation (lA.lj) let us choose I t = 0. In this case one can easily fix the null form S 
from the first equation in (IA.50 as 

s = -0 = (axx + (3x 3 ) 1Q 

Substituting this form of S into (1A.7j) we get 

dx 3 

xR x - (axx + (3x 3 )R x = -- — R. (A.ll) 

x 

Equation flA.llj) is a first order linear partial differential equation with variable 
coefficients. As we noted earlier any particular solution is sufficient to construct an 
integral of motion (along with the function S). To seek a particular solution for R one 
can make a suitable ansatz instead of looking for the general solution. We assume R to 
be of the form, 

R= {A{x) + B{x)x + C{x)x 2 Y' ^ A ' 12 ) 

where A, B and C are functions of their arguments, and r is a constant which are 
all to be determined. We demand the above form of ansatz liA.lfy) , which is very 
important to derive the Hamiltonian structure associated with the given equation, due 
to the following reason. To deduce the first integral / we assume a rational form for /, 
that is, / = gfcjl i where / and g are arbitrary functions of x and x and are independent 
of t, from which we get I x = ^ x9 ~J 9x and I x = f* 9 ~/ 9i . From (1A.5|) one can see that 
R — I x — fi3-f9± $ _ i _ f*9-f9* anc i rs = I so that the denominator of the 

9 *x J±9 J9x 

function S should be the numerator of the function R. Since the denominater of S is 
x (vide Eq. ( 1A.10I) ) we fixed the numerator of R as x. To seek a suitable function in 
the denominator initially one can consider an arbitrary form R = ^ ■ However, 
it is difficult to proceed with this choice of h. So let us assume that h(x, x) is a 
function which is polynomial in x. To begin with let us consider the case where h 
is quadratic in x, that is, h = A(x) + B(x)x + C(x)x 2 , which is a generalized version of 
the form considered in Ref. [23], where only the linear form in x was investigated (that 
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is C{x) = 0). Since R is in rational form while taking differentiation or integration the 
form of the denominator remains same but the power of the denominator decreases or 
increases by a unit order from that of the initial one. So instead of considering h to 
be of the form h = A(x) + B{x)x + C(x)x 2 , one may consider a more general form 
h = (A(x) + B(x)x + C(x)x 2 ) r , where r is a constant to be determined. The parameter 
r plays an important role, as we see below. 

Substituting (IA.12I) into (1A.11I) and solving the resultant equations, we arrive at 
the relation 



x(A x + B x x + C x x 2 ) - (axx + (3x 3 )(B + 2Cx) 



-ax(A + Bx + Cx 



c2\ 



(A.13) 



Solving equation (1A.13I) . we can fix the forms of A, B, C and r and substituting them 
into equation (1A.12I) we can get the integrating factor R. Doing so, we find 

x 



R 



(x + ^axY 
x 

(2x 2 + ax 2 x + l3x A ) 



a 1 > 8fJ 



a 2 < 



(A. 14) 



a 







where r 



4/3 



a 



± V" 2 - 8/3 



f lA.lOD and (IA.14D . respective 



. One can easily check the functions S and R given in 
y, satisfy (1A.8I) also. Finally, substituting R and S into 



the form (1A.9I) for the integral we get the integrals of motion ([2])-( 
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